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1. Introduction and Motivation 

What is the role of the center of the gauge group in the deconfinement phase transition of Yang- 
Mills theory? Answering this question is a very challenging problem and many investigations have 
been carried out in order to address this issue. The center of the gauge group has turned out to be 
relevant in two main aspects of confinement. The first one concerns the critical behaviour at the 
deconfinement phase transition and the second one is the potential mechanism of confinement in 
Yang-Mills theory. We will now examine these two features. This will allow us to discuss why it 
is interesting to study how confinement shows up in a gauge theory with a symmetry group with a 
trivial center like (7(2). 

1.1 Order of the deconfinement transition: a general approach 

In this section we discuss the order of the deconfinement phase transition in Yang-Mills theory. 
The approach is general and all possible gauge symmetry groups are taken into account ||, ||]. 

The Yang-Mills theory with gauge group G behaves differently at low and at high tempera- 
tures. The low-temperature regime is characterised by the dynamics of glueballs, which are color- 
singlet bound states of gluons. At high temperatures gluons are no longer confined into glueballs 
and form a plasma. In general, these two regimes are separated by a phase transition where the 
spontaneous breaking of a global symmetry related to the center ^(G) of G occurs. The Polyakov 
loop <J>(3c) corresponds to inserting a static quark into the gluon system and its expectation value 
(<&) = exp(— F /T) can be interpreted as a quantity related to the free energy F of a static quark at 
temperature T. The Polyakov loop is an order parameter for center symmetry breaking since it has 
the non- trivial transformation rule <J>(x)' = z $(x) under the center transformation characterised by 
the center element z £ ^(G). In the low-temperature phase, quarks are confined into hadrons and 
the free energy of a quark is infinite: (<&) = and the center symmetry is unbroken. Instead, in the 
high-temperature phase, the free energy F has a finite value: ($) ^ and the center symmetry is 
spontaneously broken. 

Some time ago Svetitsky and Yaffe [Q] have conjectured a connection between the critical 
behaviour of a gauge theory at the deconfinement transition and the critical behaviour of a scalar 
theory with a symmetry corresponding to the center of the gauge group. Integrating out the spatial 
degrees of freedom of the (d + 1) -dimensional Yang-Mills theory, one can write down a local 
effective action for <I>. This effective action describes a scalar field theory in d dimensions with a 
global symmetry ^(G). The confined phase of the Yang-Mills theory corresponds to the disordered 
phase of the scalar model since the center symmetry is unbroken. The deconfined phase has, 
instead, its counterpart in the ordered phase. 

The effective action has, in general, a complicated form and it depends on a number of param- 
eters that need to be tuned in order to match with the Yang-Mills theory and reproduce its critical 
behaviour However, if the deconfinement phase transition is second order, approaching criti- 
cally the correlation length diverges and the critical behaviour is universal. Hence, the details of 
the local effective action become irrelevant: only the center symmetry ^(G) and the dimensionality 
d of space determine the universality class. 

Svetitsky and Yaffe's conjecture assigns the center an important role in the critical dynamics 
of the Yang-Mills theory at the deconfinement phase transition. Many numerical simulations have 
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investigated the deconfinement phase transition in Yang-Mills theory and they have successfully 
checked the validity of that conjecture in the cases of second order deconfinement phase transition. 
These studies have mainly considered the Yang -Mills theory with gauge group SU (N). However, if 
one wants to investigate the relevance of the center in determining the order of the deconfinement 
phase transition, then the SU (N) groups are not the only and not even the best choice. In fact, 
since the center of SU(N) is Z(Af), when changing the group the center changes as well. This 
implies that the symmetry of the effective scalar field theory is modified and the eventually available 
universality class changes too. Hence one should widen the perspective and consider other possible 
gauge groups. In ( |1.1[ ) we list all possible gauge simple Lie groups with their center subgroups. 

( Z(2); A' odd 

Y(SU(N)) = Z(N); V(Sp(N)) = Z(2); V(SO(N))-^(Spin(N)) =1 Z(2) 2 ; N = 4k 

Z(4); N = 4k + 2 

^(G(2)) = ^(F(4))=<r(£(8)) = {ll}; ^(£(6)) = Z(3); tf(E(7)) = Z(2) 

(1.1) 

Besides the SU (N) branch, there are two other ones. The first one is the Sp(N) branch that is 
characterised by the property that all its members have the same center Z(2). The second branch 
is that of the SO(N) groups. Contrary to the other ones, these groups are not simply connected 
and their universal covering groups are Spin(N): they have center Z(2) for N odd and center 
Z(2) x Z(2) or Z(4) when N is an even or an odd multiple of 2, respectively. In this discussion 
of the deconfinement phase transition with a general group, we will assume that different groups 
with the same algebra have the same continuum limit. In particular, we will consider the universal 
covering group associated to a given algebra. Several numerical investigations concerning this 
issue and supporting this viewpoint have been carried out [^]-[13]. Finally we have 5 other groups 
which are exceptional: 3 of them, G(2), F(4) and E(S) have a trivial center while the other two, 
£"(6) and £"(7), have, respectively, the center Z(3) and Z(2). 

The groups SU (N) and SO(N) leave invariant the scalar product of ,/V dimensional complex and 
real vectors, respectively. The groups Sp(N) leave the scalar product of ,/V dimensional quaternions 
invariant. While the 5 exceptional groups leave certain forms involving octonions invariant. 

Looking at the Sp(N) branch turns out to be particularly interesting. In fact, contrary 
to SU (N), one can disentangle the size and the center of the group: one can increase the size of 
the group keeping the center Z(2) fixed. Thus, all Sp(N) Yang-Mills theories in (3+1) dimensions 
may have a continuous deconfinement phase transition in the universality class of the 3-dimensional 
Ising model. Interestingly, Sp(N) groups are another generalisation of SU (2) = Sp(l) than SU (N). 
Numerical simulations of Sp(2) and Sp(3) Yang-Mills theories in (3+1) dimensions show evidence 
for first order deconfinement phase transitions [|j]. Based also on the results available in the liter- 
ature for SU (N) groups, we conclude that the center does not play a role in determining the order 
of the deconfinement phase transition. In other words, the order of the deconfinement phase transi- 
tion is a strongly dynamical issue that can not be decided simply by symmetry arguments. On the 
other hand, as Svetitsky and Yaffe have pointed out, if the transition is second order then the center 
determines the universality class for the phase transition. 

We have conjectured []l|, ||, ||] that the size of the group determines the order of the deconfine- 
ment phase transition. At low temperature, we have the confined phase where the relevant degrees 
of freedom are glueballs, whose number is essentially independent on the size of the group. At 
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high temperature, we have, instead, a gluon plasma where the number of the relevant degrees of 
freedom increases proportionally to the size of the group. If the mismatch of the number of degrees 
of freedom between the two phases is large, then we can not move continuously from one phase 
to the other and a first order deconfmement phase transition shows up. The natural question that 
now arises concerns the finite temperature behaviour of Yang-Mills theories for which the gauge 



group has a trivial center, like e.g. G(2) Q14J]. In fact, in this case, there is no center symmetry that 
can break down and a crossover could simply connect the low- and the high-temperature regimes. 
However, according to our conjecture about the relevance of the size of the group in determining 
the finite temperature behaviour of Yang-Mills theory, we expect a first order deconfmement phase 
transition. In fact, Sp(2) has 10 generators and Sp{2) Yang-Mills theory has a first order decon- 
fmement transition in (3+1) dimensions. The group G{2) has 14 generators and so we expect G(2) 
Yang -Mills theory to have a first order transition||3|, |l5| ]. 

1.2 The center and mechanisms of confinement 

In this section we discuss the role of the center in the possible mechanism of confinement in 
non-Abelian gauge theories. 

Understanding the mechanism of color confinement stands out as one of the most challenging 
open theoretical problems of QCD. Despite the long familiarity with this phenomenon and the 
huge amount of phenomenological results, the mechanism that confines quarks inside hadrons still 
lacks a full and satisfactory explanation. The most fruitful approach is an effective description 
of the phenomenon dating back to the late '70s - early '80s JT^]-[[T9|], when it was proposed that 
confinement could be related to topological objects which, after condensation, disorder the system. 
The candidate topological objects are instantons, merons, Abelian monopoles and, in particular, 
't Hooft flux vortices, which are believed to be the most fundamental ones due to their simple 
structure. Many Monte Carlo simulations have been carried out in order to study these topological 



objects [|20|]-[]36|]. If we consider a Yang-Mills theory with gauge group G, 't Hooft flux vortices 



are present if the following inequality holds 

m{G/V(G)) ^ {1} (1.2) 

Most of the numerical investigations involve a partial gauge fixing step to detect topological 
objects. However, although the gauge-fixed approach gives interesting and intriguing results, it 
is important to emphasise that it is not on a completely solid theoretical ground due to its gauge 
dependence and the related non-locality. It should be considered as an approximate method to 
extract information about the effective mechanism of confinement. Hence, alternative approaches 
and different methods to address the problem of confinement should also be looked for. 

For example, the temperature dependence of the 't Hooft flux vortex free energy has been 
investigated. This is a gauge-invariant quantity and the numerical results confirm the relevance of 
vortices in the phenomenon of confinement. Another interesting way to investigate their role is to 
consider gauge theories without 't Hooft flux vortices and study how confinement shows up. 

The simplest pure gauge theory with only trivial 't Hooft flux vortices has the exceptional 
group G{2) as gauge symmetry Jl4|]. Numerical studies in this gauge theory can add further insight 
in the mechanism of confinement. An additional bonus that G{2) provides us is that it has SU (3) as 
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a subgroup. This allows us to move back and forth between G(2) and SU (3) Yang-Mills theories 
by exploiting the Higgs mechanism. In this way we can study how confinement in (7(2) Yang-Mills 
theory is related to the more familiar case of SU (3) Yang-Mills theory. 

2. G(2) generalities 

In this section we present the main features and properties of the exceptional group G(2). It is 
a subgroup of the real group 50(7) which has 21 generators and rank 3. The 7x7 real matrices Q. 
of 50(7) are defined by the following two properties 



det^ = l; S a b = 8 a >b> Q. aa ' &bv 

In addition to these two conditions, the matrices of G(2) satisfy also the constraint 



Tabc 



la'b'c 1 



&aa' Q-bb' Q 



(2.1) 



(2.2) 



where T a b c is a completely antisymmetric cubic tensor whose non- vanishing elements - up to index 
permutations - are 

7^127 = T\54 = 7i63 = 7235 = 7264 = 7374 = 7576 = 1 (2.3) 



The additional defining equation ( |2.2[ ) sets 7 constraints on the 50(7) matrices reducing the initial 
21 degrees of freedom to 14, which is the dimension of (7(2). Furthermore the fundamental rep- 
resentation of G(2) is 7 dimensional and so (7(2) inherits from 50(7) the property of having real 
representations. Hence G(2)-"quarks" and (7(2)-"antiquarks" live in equivalent representations and 
can not be distinguished. The rank of (7(2) is 2 and so - similarly to SU (3) - the diagrams of its 
representations can be drawn on a plane. In figure |T]a we show the fundamental representation {7}. 



1/(3 

1/2 3 


-1/2(3 
-1/(3 



(a) 




-13/2- 
1 A3- 



-1/2(3- 
-i 3 / 2- 



(b) 




-1/2 1/2 -1 -1/2 1/2 1 

Figure 1: Diagrams of the fundamental {7} (a) and of the adjoint {14} (b) representations of G(2). 



The group SU(3) is a subgroup of (7(2). When restricting to SU(3), the fundamental repre- 
sentation {7} of (7(2) becomes reducible and is given by the sum {3} © {3} © {1}. Note that this 
sum describes a real representation of SU (3) since (7(2) is a real group. This embedding can be 
made more explicit in the following basis, where 8 of the 14 (7(2) generators can be written as 

( K 

| . (2.4) 

\o o -K 



X, 



1 

7! 
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Here X a , with a = 1,2, . . . , 8, are the usual 3x3 Gell-Mann matrices. The SU (3) embedding can 
also be read off directly from the diagram in figure [l]a, where one can easily identify the triangles 
of the {3}, the {3} and the trivial representations of SU (3). 

The fact that SU (3) is a subgroup of G(2) can be used to derive an explicit representation 
of the matrices [p7|]. In fact, a (7(2) matrix, O, can be expressed as the product of two matrices 
and i2°: the former, % ', belongs to the SU (3) subgroup while the latter, is in the coset 
G(2)/SU(3) ~ S 6 . The matrices and ^ are given by 



/ C(K) n(K)K 



-li{K)tf L_m —^i[K)K T I ; flf(ff) 
\ D(K) n(K)K* 



(u 
1 

u* 

(2.5) 

where K is a 1 x 3 complex vector and U is a 3 x 3 SU (3) matrix. The numbers x{K) and }l(K) are 
given by x(K) = \ \K\ | 2 and }l{K) = \/2/(l +x{K)), while the 3 x 3 matrices C{K) and D(K) are 




V 



C(K) = - < 11 ^ ; D(#) = T where , ^ 

v ; A\ A(l+A)/' v 7 A A 2 = £ a pyK r ; A = y'T+^C^T 

(2.6) 

Using this parametrisation one can see that £2 depends on 14 real parameters since the vector K and 
the SU (3) matrix £/ depend, respectively, on 6 and 8 real parameters. 

The group (7(2) has 14 generators and so the adjoint representation is 14 dimensional: its 
diagram is shown in figure |l|b. Again, restricting to the SU (3) subgroup, the adjoint representation 
becomes reducible and it is given by the sum {8} © {3} © {3}. Thus, from the viewpoint of the 
SU (3) color degrees of freedom, 8 of the G(2)-"gluons" behave as the 8 gluons and the remaining 
6 split into a triplet and an antitriplet behaving like vector quarks and antiquarks. Similar to the 
case of the fundamental representation, the reduction of the adjoint representation of (7(2) down 
to SU (3) can also be read off from the diagram in figure [T]b. Indeed we identify the two triangles 
corresponding to the {3} and the {3} and the hexagon of the adjoint representation {8} of SU (3). 

The fact that both Gil) and SU(3) have rank 2 and that SU(3) C G(2), implies that the center 
of (7(2) is a subgroup of Z(3), the center of SU (3). The property of having real representations 
then makes the center of (7(2) trivial. This feature - together with the property that (7(2) is its own 
universal covering group - leads to the absence of non-trivial 't Hooft sectors: 

m (G(2)/tf(G(2))) = 11 (2.7) 

The triviality of the center also implies that all G(2) representations mix together in the tensor 
product decompositions and there is no property similar to triality that splits the SU (3) represen- 
tations into 3 branches. Hence, in contrast to SU(3), a heavy G(2)-"quark" can be screened by 
G(2)-"gluons" and the string breaks already in the pure glue theory without the need for dynamical 
G(2)-"quarks". In particular, from the tensor product decomposition 

{7} ©{14} ©{14} ©{14} ={!} + ••• (2-8) 
one concludes that 3 G(2)-"gluons" are sufficient to screen a G(2)-"quark". 
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Finally, as a last feature, we consider the following homopoty groups that tell us which topo- 
logical objects can be expected in G(2) Yang-Mills theory 

7T 3 (G(2)) = Z; tt 2 (G(2)/£/(1) 2 ) = ZxZ (2.9) 

Hence, there are two types of monopoles just like in the SU (3) Yang-Mills case and we have 
instanton solutions leading eventually to non-trivial 0-vacuum effects. 



3. G(2) Yang-Mills theory 

In this section we discuss the Yang -Mills theory with gauge symmetry group G(2), presenting 
results of numerical simulations at finite temperature. 

As we have seen in the previous section, G(2) Yang-Mills theory describes the interaction of 
14 G (2) -"gluons". As any non-Abelian pure gauge theory in 4 dimensions, it is asymptotically free 
with a non-perturbatively generated mass-gap. Thus we expect a confining phase at low tempera- 
tures but with a vanishing string tension since the string between two static G(2)-"quarks" breaks 



by production of dynamical G(2)-"gluons". Equation ( |2.8| ) shows that at least 3 G(2)-"gluons" are 
needed to screen a G(2)-"quark". As we have seen in the previous section, restricting to the SU (3) 
subgroup, 8 of the 14 G(2)-"gluons" are related among themselves as the usual gluons of QCD and 
the remaining 6 split up into a triplet and an antitriplet. Hence, just like quarks in QCD, these last 
6 G(2)-"gluons" explicitly break the center symmetry Z(3) of SU (3). Thus, in contrast to SU (3) 
Yang-Mills theory, in G(2) Yang-Mills theory we expect the phenomenon of confinement to resem- 
ble more closely the one of QCD but yet without the complications related to dynamical fermions. 

Since G(2)-"gluons" screen G(2)-"quarks", the Wilson loop always has - up to short dis- 
tance effects - a perimeter law and its behaviour can not be considered as an order parameter for 
confinement. In other words, the string tension, defined as the ultimate slope of the static G(2) 
"quark"-"quark" potential, vanishes due to string breaking and it can not be used to characterise 
the low-temperature phase of G(2) Yang-Mills theory. One can consider the Fredenhagen-Marcu 
order parameter [^] which is a quantity probing whether the theory has a long range interaction or 
not, i.e. if we are in a massless Coulomb phase or in a massive confining/Higgs phase. By a strong 
coupling expansion one can show that the theory is indeed in a massive phase. 

At finite temperature, due to the triviality of the center of G(2), we expect a different behaviour 
than in SU (3) Yang-Mills theory, where confined and deconnned phases differ by the way the cen- 
ter symmetry is realized. Since in G(2) Yang-Mills theory no global center symmetry can break 
down, we can almost exclude a second order deconfinement phase transition and a crossover can 
simply connect the low and high temperatures. However, also a first order transition at finite tem- 
perature can take place, separating the low- and high-temperature regimes. Although a crossover 
may seem the most natural possibility, according to the conjecture we have discussed in section 
|1.1| , we expect a first order transition. In fact, in the low-temperature regime, the dynamics is well 
described in terms of color-singlet glueball states, whose number depends mildly on the size of the 
group, i.e. on the number of gluons. On the other hand, at high temperatures, the large number 
of gluons, 14, becomes very relevant since we have a gluon plasma. Hence the large mismatch in 
the number of the relevant effective degrees of freedom between low and high temperatures, sug- 
gests a discontinuous behaviour at finite temperature with a first order phase transition. However 



017/7 



Confinement and the center of the gauge group 



Michele Pepe 



the presence of a finite temperature phase transition and its order are dynamical features that can 
be properly addressed only by numerical simulations. In order to investigate these issues we have 
performed Monte Carlo simulations of G(2) Yang -Mills theory in 4 dimensions on the lattice. 

The G(2) Yang-Mills theory on the lattice is constructed in the usual way in terms of link 
matrices, U x> ^, which are elements of the group G(2) in the fundamental representation {7} . We 
consider the Wilson action, which is given by 



S YM [U] 



Tr 



x,[i<v 



U x uU x 



(3.1) 



where g is the bare gauge coupling. As a first step, we have measured the action density at many 
values of the gauge coupling in order to check for the presence of a possible bulk phase transition 
separating the strong coupling from the weak coupling regime, where the continuum limit can be 
approached. In figure [2|a we show the numerical results. The red and green points correspond to 
a fine-grained hysteresis cycle in the region where a bulk transition could occur. This approach 
has been used in order to enhance the efficiency in detecting coexisting phases in case of a bulk 
transition. Our results indicate that no bulk transition separates the strong and the weak coupling 
regimes. The red and the blue lines are, respectively, the analytic perturbative expansions at strong 
and weak coupling. 



TtUnp 




4> 




Vis 1 ) 



Figure 2: (a): Action density Tr t/rj/7 as a function of the gauge coupling 1 /(g 2 ). (b): Polyakov loop Monte 
Carlo history from a numerical simulation on a 20 3 x 6 lattice at 1 /(g 2 ) = 9.765. 

Since no bulk phase transition interfering with the deconfinement phase transition has been 
found, we have looked for the presence of a transition at finite temperature. Consistent with our 
expectation and contrary to the argument based on the absence of a center symmetry, we have found 
that the G(2) Yang-Mills theory has a deconfinement phase transition at finite temperature 
In figure ||b we show the Monte Carlo history of the Polyakov loop in the critical region. A 
first order deconfinement phase transition can be clearly observed, with many tunneling events 
between two coexisting phases. It is important to note that, since no symmetry gets broken, no 
order parameter can be defined. However discontinuities in physical quantities, for instance in the 
free energy of a static G(2)-"quark" or in the specific heat, can be used to unambiguously detect this 
finite temperature phase transition. A detailed analysis requires a finite size scaling study and the 
continuum limit extrapolation. The value of the lattice temporal extent, N t = 6, and the set of gauge 
couplings considered in this study have been chosen in order to be close to the continuum limit. 
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Figure 3: Polyakov loop probability distributions in the critical region of the deconfinement phase transition 
in the (3+l)-dimensional G(2) Yang-Mills theory. The temperature increases from left to right. 



In figure || we show the probability distributions of the Polyakov loop in the critical region: 
the temperature increases from left to right and the switching from one phase to the other can be 
clearly observed. Notice that the Polyakov loop in the low-temperature, confined phase is peaked 
around zero. However, due to string breaking, we know that the free energy of a static G{2)- 
"quark" is always finite. Hence, the very small value we observe here - hardly distinguishable 
from zero - is related to the fact that the energy necessary to pop out the G(2)-"gluons" that screen 
the G(2)-"quark" from the vacuum is quite large. 



4. The high-temperature effective potential 

In the previous section we have presented the results of numerical simulations in G(2) Yang- 
Mills theory, showing the presence of a finite temperature, first order deconfinement phase transi- 
tion. In order to have a better understanding of the deconfined phase, we discuss here the analytic 
perturbative computation of the effective potential for the Polyakov loop at high temperature. 

Perturbative computations at zero and finite temperature differ in an important aspect. At zero 
temperature, the temporal component of the gauge field can be eliminated by a gauge transfor- 
mation defined up to static gauge transformations. At finite temperature, instead, the temporal 
component can not, in general, be gauged to zero but only to a static value due to the compactified 
temporal direction. Hence, perturbative expansions around background gauge fields with different 



static temporal components are physically different at finite temperature [ 39j , [4( 

Here we compute the one-loop free energy in the continuum. We perturbe around a general 
static background gauge field given by 

A* = 8^{V2 d 1 A 3 + V6d 2 Az) (4.1) 

where dy and 62 are the two phases characterising an Abelian G(2) matrix. Note that evaluating the 
effective potential for the background field A^ is equivalent to computing the effective potential for 
the Polyakov loop since it is simply the time-ordered integral of A^ along the temporal direction. 
We now introduce a gauge fixing condition defined by D^A^ = 0, where the covariant derivative is 
given by D?f = d^5 ab + gf abc A c ^, with / being the G(2) structure constant. After integrating in 
the ghost field %, the continuum Lagrangian takes the form 

2 = \g%v + ^(Df l A fl ) 2 + xDf l D^x (4-2) 
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where a is a gauge fixing parameter and the field strength is defined as usual G"^ v = d^A" — 
d v A'^ — gf abc A h ^A b v . We now decompose the gauge field in the sum of the background field and of 
the quantum fluctuations : A^ = + A^ . Retaining in the Lagrangian only terms that are at 



most quadratic in A^ , the partition function can be computed and the free energy W is then given by 

(8, V (-D B ) 2 ) + (1 - ^K»v] } " Tr [log(-(D B ) 2 )] (4.3) 

Using some formulas [ fill ] to manipulate the previous expression, the one-loop free energy can be 
explicitly calculated and the result is 



w[e l ,e 2 } = -Tr{io g 



4n 2 T 4 



30 



Q(d u d 2 ) 



2k 



(4.4) 



where Bn{x) = — l/30+x 2 (x— l) 2 is the 4 th Bernoulli polynomial with the argument defined mod- 
ulo 1, l/T is the length of the compactified temporal direction and d\ = g6\/T and d 2 = gd 2 /T . 
The C/'s are given by 



Ci(0 1 ,0 2 )=20 1 ; C 2 (0 1) 2 ) = 1 +30 2 ; 
Ca{0\,Q 2 ) =2B 2 ; C 5 (di,d 2 ) = di — d 2 ; 



c~3(Qi,Q 2 ) — Qi -36) 2 ; 
C6(di,d 2 ) = di + 6 2 ; 



(4.5) 



Close to the trivial background field (6> l5 G 2 ) = (0,0), we have W = — (14/45)7T 2 r 4 , which corre- 
sponds to an ideal gas of G(2)-"gluons"; note the number 14 in the numerator corresponding to 
the number of G(2)-"gluons". Finally, we point out that, since SU (3) is a subgroup of G(2) with 
the same rank, the high-temperature effective potential in SU (3) Yang-Mills theory [|{], ^TJ can be 
immediately obtained from ( |44| ) performing the sum only up to 3. In figures |]a and^|b we show 
the contour plots of the one-loop effective potential for the Polyakov loop at high temperature in 
SU (3) and (7(2) Yang-Mills theories, respectively. 

(b) 

2Kv ;im^^»|o^ 



0, o 




2 o 



-3-2-10 1 2 3 




0i 

Figure 4: Contour plots of the one-loop high-temperature effective potential for the Polyakov loop for SU (3) 
(a) and G(2) (b) Yang-Mills theories. 



The formula (1.4) is, somehow, a counter-intuitive result. In fact, despite the fact that the 
center of G(2) is trivial, the center Z(3) of the SU (3) subgroup plays a role in the effective po- 
tential of G{2) at high temperature (remember that c to{G(2)) C c if(SU(3))). Indeed we have an 
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absolute minimum corresponding to the trivial center element (61,62) = (0,0) but we also have 
local minima corresponding to the non-trivial Z(3) center elements (61,62) = (0,±27i/3) and 
(61,62) = (±7T,±7r/3). The non-trivial Z(3) center elements are metas table minima since they 
are separated from the absolute minimum by a free energy gap that grows proportionally to the 
volume: hence they are not relevant in the thermodynamic limit and in determining the bulk be- 
haviour of the (7(2) Yang-Mills theory. However, numerical simulations are carried out on lattices 
of finite size and so there is a finite probability of tunneling to one of these metastable minima. 
Then, depending on the algorithm used to update the system, these metastable states can have long 
lifetimes and significantly bias the numerical results. 

5. G(2) Yang-Mills theory with a Higgs field 

In this section we discuss the connection between confinement in G(2) Yang-Mills theory and 
in the familiar case of SU (3) Yang-Mills theory. In fact, by the Higgs mechanism, we can break the 
G(2) gauge symmetry down to SU (3) and study how the well-known form of confinement in SU (3) 
Yang-Mills theory reemerges as the breaking of the G(2) symmetry becomes stronger and stronger. 

The gauge symmetry breaking of G(2) Yang-Mills theory down to SU (3) can be accomplished 
by adding a Higgs field in the fundamental representation {7} of G(2). Restricting to the SU (3) 
subgroup, the adjoint representation {14} is reducible and splits into the sum {8} © {3} © {3}. 
Hence, the 8 G(2)-"gluons" which are related among themselves like the gluons, stay massless 
while the other 6 pick up a mass Mq proportional to the expectation value v of the Higgs field. 
Changing v, we can then change the mass of the 6 massive G(2)-"gluons". If Mq is not too large 
compared to Aqqd, they participate in the dynamics. However, as Mq becomes larger and larger, 
they progressively decouple and, at the end, they are removed from the dynamics and we are 
left with the SU (3) Yang-Mills theory. Thus a Higgs field in the fundamental representation {7} 
provides us with a handle we can use to interpolate between G(2) and SU (3) Yang-Mills theories. 

In SU (3) Yang-Mills theory the behaviour of the Wilson loop is an order parameter for con- 
finement: the string tension is non-zero in the confined phase and vanishes in the deconnned phase. 
On the contrary, in G(2) Yang-Mills theory, the string tension is always zero. We now discuss how 
the Higgs mechanism relates these two scenarios. Let us first consider the pure glue case. As we 
have discussed in section [| the breaking of the string between two static G(2)-"quarks" happens 
due to the production of 3 pairs of G(2)-"gluons". Hence, the string breaking scale is related to the 
dynamical mass of the 6 G(2)-"gluons" popping out of the vacuum. When we switch on the inter- 
action with the Higgs field, 6 of the G(2)-"gluons" start to pick up a mass for the Higgs mechanism. 
Note now that the two G(2)-"quark"/"gluon" bound states resulting from string breaking, must be 
both G(2)-singlets and SU (3)-singlets. Then, since the fundamental representation {7} reduces 
to the sum {3} © {3} © {1} when restricting to the SU (3) subgroup, it follows that, among the 6 
G(2)-"gluons" breaking the string, there must be a pair of the massive ones. Thus the string break- 
ing scale depends on Mq and the larger Mq goes, the larger is the distance where string breaking 
occurs. The picture of the unbreakable SU (3) string then reemerges. When the expectation value 
of the Higgs field is sent to infinity, so that the 6 massive G(2)-"gluons" are completely removed 
from the dynamics, also the string breaking scale is infinite. Hence, we fully recover the familiar 
SU (3) linearly rising, confining potential with a non-vanishing value for the string tension. 
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In order to have a more quantitative understanding of the relation between confinement in 
G(2) and SU (3) Yang-Mills theories, we have performed numerical simulations in G(2) Yang- 
Mills theory with a Higgs field in the fundamental representation {7}. In particular, the interesting 
aspect we want to address is the connection between the finite temperature critical points in the two 
pure gauge theories. We consider a Higgs field (p of unit length; the lattice action looks as follows 

Shym [U, <p] = S YM [U]-kY, <PxUx,h 9x+fi (5- 1) 

X,jl 

where k is the parameter characterising the interaction between the Higgs field and the G(2)- 
"gluons". In our numerical study we have considered the same temporal extent, N t = 6, that we 
have used in the Monte Carlo simulations of G(2) Yang-Mills theory. The aim of our numerical 
simulations has been to understand the phase diagram in the parameter space (l/(g 2 ),fc). The 
result of this study is summarised in figure [|. 



K 



\ 




c 


1 


1. 



Figure 5: Phase diagram in the parameter space (1/ (g 2 ), k). The K — and K — °° axes correspond to the 
G(2) and SU (3) Yang-Mills theories respectively. The l/(g 2 ) = °° limit is the SOil) spin model. 

Let us first discuss some particular points in the phase diagram. For K = 0, the Higgs field 
decouples and we have the pure G(2) Yang-Mills theory: the point A in figure || is the correspond- 
ing critical coupling of the deconfinement transition. For fc = oo, due to the Higgs mechanism, the 
6 massive G(2)-"gluons" are removed from the dynamics and the G(2) gauge symmetry is com- 
pletely broken down to SU (3). The point B in the phase diagram is the critical coupling of the 
SU (3) Yang-Mills theory at N t = 6 (note that a rescaling factor 7/6 has been taken into account 
for the different normalisation of the coupling). For l/(g 2 ) = °°, the gauge degrees of freedom 
are completely frozen and we have a spin model with global 50(7) symmetry. This model has an 
order/disorder phase transition that is denoted by the point C in the phase diagram. 

Starting from A, B and C, we have carried out numerical simulations in order to investigate 
how these points are related. In particular, our main interest is to study the relation between A 
and B. However let us first discuss the point C. We have not found any indication of a critical line 
ending at the point C. Nevertheless, since it was not our main interest studying this part of the phase 
diagram, we can not rule out a second order line or a weak first order one that could, eventually, 
show up considering larger lattices. For this reason in figure || we have plotted a dotted line ending 
at the point C. Let us now consider the point A. As we enter the phase diagram switching on the 
gauge-Higgs interaction, we observe that the first order transition is still present and that the critical 
coupling slowly shifts towards smaller values. In figure ^a we show the probability distribution for 
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the Polyakov loop at 3 values of the gauge coupling - increasing values from bottom to top - at 
fixed gauge-Higgs coupling K = 1.3. Two coexisting phases show up and their relative relevance 
switches as the temperature increases: this gives a clear signal of a first order transition. 

(a) ^ (b) ^ (c) 



- 
0.(12 - 




Figure 6: Polyakov loop probability distributions at K = 1.3 (a), K = 1.5 (b) and K = 4.0 (c). The gauge 
coupling 1 / (g 2 ) increases from bottom to top. 

However, as we consider larger values of K, the first order transition weakens and then the 
phase transition is washed out to a crossover. For instance, in figure |6|b, we show the probability 
distributions of the Polyakov loop at 3 different values of the gauge coupling (increasing values 
from bottom to top) at K = 1.5. The distributions are quite broad and we do not find any indication 
of a phase transition. However, we can not exclude that, performing numerical simulations on 
much larger volumes, a weak first order transition could show up. 

If we further increase the parameter K, the phase transition starts to show up again and it 
stays there up to fc = oo, where the limit of the SU (3) Yang-Mills theory is attained. Note that the 
probability distribution of the real part of the Polyakov loop at the deconfinement phase transition 
in SU (3) Yang-Mills theory is characterised by 3 peaks: one corresponding to the confined phase 
and the other two associated with the deconfined one. One of the peaks of the deconfined phase 
is related to the trivial Z(3) center element while the other one is the projection on the real axis 
of the two non- trivial Z(3) center elements. We expect that a similar 3 -peak structure should also 
progressively show up in the G(2) gauge-Higgs system when the parameter K is large enough. 
However, it is important to point out that, as long as K has finite values, the peak corresponding to 
the non-trivial Z(3) center elements is metastable and there is a free energy difference proportional 
to the volume with respect to the other deconfined peak. Hence, the second deconfined peak plays 
no role in the critical behaviour since, in the thermodynamic limit, it would be absent. Nevertheless, 
in a finite volume, as it is the case when performing numerical simulations, the second deconfined, 
metastable peak should be observable for sufficiently large K. Indeed, as we can see in figure^, the 
3 peak structure shows up. Similar to the two previous plots, we report the probability distribution 
for the Polyakov loop at 3 different values of the gauge coupling, at fixed K = 4.0. 
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6. Conclusions 

Understanding the phenomenon of confinement and the features of the deconfinement phase 
transition is an open, challenging problem. In this paper we have discussed the role played by the 
center of the gauge group: we have approached the problem from the most general perspective 
of considering all possible choices for the gauge symmetry group. Then we have focused our 
attention on the case of G(2). The group G(2) is the smallest, simply-connected group with a 
trivial center and G{2) Yang-Mills theory is the simplest pure gauge theory with no non-trivial 
't Hooft flux vortices. We have performed Monte Carlo simulations in G{2) Yang-Mills theory in 
(3+1) dimensions, finding numerical evidence for a first order deconfinement phase transition. This 
result is consistent with our conjecture that the size of the gauge group and not the center plays a 
relevant role in determining the order of the deconfinement transition. 

By exploiting the Higgs mechanism, one can study the connection between the deconfinement 
phase transition of G{2) Yang-Mills theory and the more familiar one of SU (3) Yang-Mills theory. 
The numerical results indicate that the first order deconfinement transition of G{2) Yang-Mills 
theory weakens as the interaction with the Higgs field is switched on. This supports our conjecture 
on the relation between the size of the group and the order of the deconfinement transition. In fact 
the Higgs mechanism removes progressively a number of degrees of freedom from the dynamics. 

As the gauge-Higgs interaction is furtherly increased, a first order transition appears again and 
it stays there until the SU (3) limit is reached. We would interpret this phase transition as an effect 
of the lack of a universality class with aZ(3) symmetry in 3 dimensions and not as a results of the 
mismatch of the number of the relevant degrees of freedom in SU (3) Yang-Mills theory between 
the confined and the deconfined phases. 
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